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t— -} , Abstract 

o; 

(y-) i For every list of integers X\, . . . , x m there is some j such that X\ + ---+Xj — £j+i — ■ • • — x m « 0. 

So the list can be nearly balanced and for this we only need one alternation between addition 
and subtraction. But what if the Xi are fc-dimensional integer vectors? Using results from 

Q topological degree theory we show that balancing is still possible, now with k alternations. 

This result is useful in multi-objective optimization, as it allows a polynomial-time corn- 
ed putable balance of two alternatives with conflicting costs. The application to two multi-objective 
i 1 1 optimization problems yields the following results: 

• A randomized 1 /2-approximation for multi-objective maximum asymmetric traveling sales- 
^. ■ man, which improves and simplifies the best known approximation for this problem. 

• A deterministic 1 /2-approximation for multi-objective maximum weighted satisfiability. 

in ; 

C-^ ■ 1 Introduction 

o: 
o: 

Balancing Sums of Vectors. Suppose we are given a sequence of goods gi, ■ ■ ■ ,g m , each of which 
has a value, a weight, and a size. Is it possible to distribute the goods on two trucks such that 
the loads are nearly the same with respect to value, weight, and size? We show that this is always 
possible by a very easy partition: For suitable indices i,j, k, I, assign gi,g, L+ i, . . . ,gj, gk,gk+i, ■ ■ ■ ,gi 
to the first truck and the remaining goods to the second one. In general, if the goods have 2k 
criteria (value, weight, size, . . . ), then there exist k intervals of goods such that the goods inside 
and the goods outside of the intervals are nearly equivalent with respect to all criteria. 

More formally, let x\, . . . ,x m £ N 2fc be vectors of natural numbers that represent the criteria of 
each good, and let z G N be an upper bound for these vectors (i.e., Xi < z for all %). Lemma 12.61 
provides intervals I\, . . . , !& C N such that for / = |J i=1 Ii, 

—4kz < 2_, x i ~ /, x i ^ 4A;z, 

where the < hold with respect to each component. The same is true if x\, . . . , x m € Z are vectors 
of integers, where — z < x% < z for all i (Corollary I2.T|) . The proofs of these balancing results are 
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based on the Odd Mapping Theorem, a result from topological degree theory, which we apply in a 
discrete setting. The discretization is responsible for the term Akz, which is caused by a rounding 
error that unavoidably occurs at the boundaries of the intervals Ji, . . . , ifc. 

The simplicity of the desired partition (i.e., a union of k intervals) is important for the application 
of our balancing results. Algorithmically, it means that for fixed dimension 2k, the right choice for 
the intervals I±, . . . , Ik can be found by exhaustive search in time polynomial in to. 



Multi-Objective Optimization. Many real-life optimization problems have multiple objectives 
that cannot be easily combined into a single value. Thus, one is interested in solutions that are 
good with respect to all objectives at the same time. For conflicting objectives we cannot hope 
for a single optimal solution, but there will be trade-offs. The Pareto set captures the notion of 
optimality in this setting. It consists of all solutions that are optimal in the sense that there is 
no solution that is at least as good in all objectives and better in at least one objective. So the 
Pareto set contains all optimal decisions for a given situation. For a general introduction to multi- 
objective optimization we refer to the survey by Ehrgott and Gandibleux [EGOQ] and the textbook 
by Ehrgott [Ehr05j . 

For many problems, the Pareto set has exponential size and hence cannot be computed in poly- 
nomial time. Regarding the approximability of Pareto sets, Papadimitriou and Yannakakis [PYOOj 
show that every Pareto set has a (1 — ^-approximation of size polynomial in the size of the instance 
and 1 /e (for the formal definition of approximation see section [3~Tj) . Hence, even though a Pareto set 
might be an exponentially large object, there always exists a polynomial-size approximation. This 
clears the way for a general investigation of the approximability of Pareto sets of multi-objective 
optimization problems. 

In general, inapproximability and hardness results directly translate from single-objective opti- 
mization problems to their multi-objective variants. On the other hand, existing approximation 
algorithms for single-objective problems can not always be used for multi-objective approxima- 
tion. Using our balancing results we demonstrate a translation of single-objective approximation 
ideas to the multi-objective case: We obtain a randomized i/^-approximation for multi-objective 
maximum asymmetric TSP and a deterministic ^-approximation for multi-objective maximum 
weighted satisfiability. 



Traveling Salesman Problem. The (single-objective) maximum asymmetric traveling salesman 
problem (MaxATSP, for short) is the optimization problem where on input of a complete directed 
graph with edge weights from N the goal is to find a Hamiltonian cycle of maximum weight. Enge- 
bretsen and Karpinski [EK01] show that MaxATSP cannot be (319/320 + e)-approximated (unless 
P=NP). In 1979, Fisher, Nemhauser and Wolsey [FNW79] gave a ^-approximation algorithm for 
MaxATSP (remove the lightest edge from each cycle of a maximum cycle cover and connect the 
remaining parts to a Hamiltonian cycle). Since then, many improvements were achieved and the 
currently best known approximation ratio of 2 /3 for MaxATSP is given by Kaplan et al. [KLSS05J. 

The fe-objective variant fc-MaxATSP is defined analogously with edge weights from N fc . The 
hardness results for MaxATSP directly translate to its multi-objective variant (just set all but one 



component of the edge weights to a constant), but algorithms have to be newly designed. Blaser et 
al. [BMP08J show that fc-MaxATSP is randomized (t+x — e)-approximable. This was improved by 
Manthey [Man09j to a randomized (^ — ^-approximation for all (fixed) numbers of criteria. Both 
algorithms extend the cycle cover idea to multiple objectives. With a surprisingly simple algorithm 
we improve the approximation ratio to 1 /2. 



Satisfiability. Given a formula in conjunctive normal form and a non-negative weight in N for 
each clause, the maximum weighted satisfiability problem (MaxSAT, for short) aims to find a 
truth assignment such that the sum of the weights of all satisfied clauses is maximal. The first 
approximation algorithm for MaxSAT is due to Johnson |Joh74| . He proved an approximation ratio 
of (2 r - 1 )/(2 r ) for formulas where each clause has at least r literals. His work showed that the general 
MaxSAT problem is 1 /2-approximable. Yannakakis [Yan94] improved the approximation ratio of 
MaxSAT to 3 /4, and Goemans and Williamson [GW94] subsequently gave a simpler algorithm with 
essentially the same approximation ratio, and later [GW95] improved the approximation ratio to 
0.758. Further improvements followed, and the currently best known approximation ratio of 0.7846 
is due to Asano and Williamson [AW02]. Regarding lower bounds, Papadimitriou and Yannakakis 
[PY91J show that MaxSAT is APX-complete. Furthermore, by Hastad |Has97| . MaxSAT cannot 
be approximated better than 7 /s, unless P=NP. 

Only little is known about the multi-objective maximum weighted satisfiability problem 
(fc-MaxSAT, for short), where each clause has a non-negative weight in N for some fixed k > 1 
and where we wish to maximize the weight of the satisfied clauses. Santana et. al. [SBLL09] apply 
genetic algorithms to a version of the problem that is equivalent to fc-MaxSAT with polynomially 
bounded weights. To our knowledge, the approximability of fc-MaxSAT has not been investigated 
so far. 

Using our balancing results, we can transfer a simple idea from single-objective optimization to 
the multi-objective world: For any truth assignment, the assignment itself or its complementary 
assignment satisfies at least one half of all clauses. We obtain a (deterministic) l /2- approximation 
for fc-MaxSAT, independent of k. 



2 Balancing Results 

2.1 Preliminaries 

Let a, b £ R. We call a function /: [a,b] — > R integrable, if it is Lebesgue-integrable on [a, b]. 
This is especially the case for bounded functions / with only finitely many points of discontinuity. 
A function g: [a,b] — > W 1 is componentwise integrable, if all projections gi are integrable and 
in this case we write J g(x) dx as abbreviation for the tuple ( f gi(x) dx, . . . , J g n (x) dx). For 
x = (xi, . . . , x n ), y = (yi, ... , y ri ) € W 1 we write x < y if X{ < yi for all i G {1,2,..., n}. For a 
set A C M. n , A denotes the (topological) closure of A, and dA denotes the boundary of A. The set 
A C M n is symmetric if x £ A ^=^> —x G A for all x G M. n . 



For bounded, open sets D C R n , continuous functions p: D — > R n and points p € R n \ p(dD) 
the integer d(tp,D,p) is called the Brouwer degree of ip and -D at the point p. We will not define 
it here, but we note that it captures how often p is "covered" by <p(D), counting "inverse" covers 
negatively, and that it generalizes the winding number in complex analysis. 



2.2 Analytical Version 

We apply the following theorems from topological degree theory to get the analytical version of our 
balancing results. 

Theorem 2.1 ( |Llo78t Theorem 2.1.1]). If D C R n is bounded and open, <p: D — > R n is continuous, 
p $l p{dD) and d(p, D,p) ^ 0, then p € p(D) . 

Theorem 2.2 (Odd Mapping Theorem, [Llo78, Theorem 3.2.6]). Let D be a bounded, open, sym- 
metric subset ofW 1 containing the origin. If <p: D —> R™ is continuous, ^ p(dD), and for all 
x € 3D it holds that -pSr 7^ i(_^i > then d(ip, D,0) is an odd number (and in particular not zero). 

Corollary 2.3. Let D be a bounded, open, symmetric subset of R n containing the origin. If 
ip: D — ^ R n is continuous and for all x S dD it holds that p(—x) = —p(x), then € <p(D). 

Proof. Assume that ^ 'p(D). From p(— x) = —p{x) for x € dD it follows that the inequal- 
ity condition of Theorem 12.21 is fulfilled (note that ^ p(dD)) and thus d(p,D,0) 7^ and by 
Theorem 12.11 € *p(D). This is a contradiction. □ 

Lemma 2.4. Let n > 1, a,b £ R ; and h: [a,b] — > R 2n be componentwise integrable. There exist n 
closed intervals I\, . . . , I n C [a, b] such that for I = I\ U • • • U I n , 



h(x) dx = / h{x) dx. 
I [a,b]\I 

Proof. Observe that it suffices to show this for [a, b] = [0, 1]. Define T = {(t\,t2, • • • , £2™) £ R 2n 
Ya=i l*il ^ !} an d for every t = (t 1 ,... ,t 2n ) G T, let 



u 


'k-1 




Kfc<2n 


-i=i 


i=l 


t k >0 







and 

/:T^R 2n , f{t)= J h{x)dx- I h(x)dx. 
h [0,1] \j t 

By the formal definition, /^ is a union of (at most) 2n closed intervals. However, it can always be 
written as a union of at most n closed intervals, by merging adjacent intervals. 

We now want to show that & f(T) by applying Corollary 12.31 to ip = f and D being the interior 
of T. D is obviously a bounded, open, and symmetric subset of R 2n containing the origin. The 
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Figure 1: Illustration of the set It for some value of t = (ti, . . . , t$), where t\, t2, t± and is are positive and 
£3 j ^5 j ^6 an d £7 are negative. 

function / is continuous because of the fundamental theorem of calculus for the Lebesgue integral 
and the fact that the endpoints of the intervals in It depend continuously on t. Furthermore, for 
any t G dD we get that there are only finitely many points in [0, 1] which are not in exactly one 
of the sets I-t and 1% and thus f(—t) = —f(t) since these finitely many points have no influence 
on the values of the integrals. Since all preconditions of the corollary are fulfilled, we get € f(T) 
and thus there exists some t E T such that 



h(x) dx = / h(x) dx. 

h [0,1] \i t 

As already noted, It can be written as a union of at most n closed intervals. We obtain a union of 
exactly n intervals by adding intervals [a, a]. □ 

Lemma 2.5. Let n>l, a,b £ R, and f,g: [a,b] — > R 2n be componentwise integrable. There exist 
n closed intervals I\, . . . , I n C [a, b] such that for I = I\ U ■ ■ ■ U I n , 

I f{x)dx+ I g(x)dx = - I f(x)+ g(x)dx. 
I [a,b]\I [a,b] 

Proof. Applying Lemma 12.41 to h(x) = f(x) — g(x) yields some / C [a, b] that is the union of n 
closed intervals in [a, b] such that 



h(x) dx = / h(x) dx 

I [a,b]\I 

f{x)-g{x)dx= J f{x)-g{x)dx 

I [a,b]\I 



/ f{x) - g(x) dx+ j g(x) - f{x) dx = 



I [a,b]\I 



2 I f(x)dx + 2 / g(x)dx= / f{x)+g(x)dx 

[a,b]\I [a,b] 

f(x)dx+ / g(x)dx = - / f(x)+g(x)dx. 

I [a,b]V [a,b] 

Note that (*) is obtained by adding J, fi , f{x) + g(x) dx to both sides. □ 
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Figure 2: Graphs of the functions / and g used in the proofs of the Lemmas 12.61 and 12.81 

2.3 Discretization of the Analytical Results 

Now we discretize the analytical results which causes a rounding error that cannot be avoided. 

Lemma 2.6. Let n, m > 1 and X\, . . . , x m ,y±, . . . , y m , z G N 2n such that Xi < z and yi < z for 
all i. There exist natural numbers 1 < a\ < b\ < 02 < 62 < • • • < o n < 6 n < m swc/t £/ia£ /or 
^ = U"=i{«i> Oj + 1, . . . , 6» - 1}, 

-. m -. m 

-2nz+ -^(xj + yj) < ^2 x i + ^2vi < 2nz + - ^(z; + y^). 



j=l 



iei 



Hi 



i=X 



Proof. For the proof it is advantageous to start the indices of Xi and yi at 0. We first define 
two functions / and g that distribute the values Xq, . . . , x m -i,yo, . . . , y m -i G N 2n equally over the 
interval [0,m), and then we apply Lemma [231 Let f,g: [0,m] — > M? n such that 



/(*) 



-"] ■-/' 



if t € [i, * + V 2 ) 



(0, . . . , 0) otherwise 



and 



0(*) 



2i/i 



if t € [i + l /2, i + 1) 



(0, . . . , 0) otherwise. 



Figure [2] shows the graph of / and g. Note that both functions are componentwise integrable. 
Moreover, for i £ {0, . . . , m — 1}, 

i-i+l fi+l 

(1) 



/ /(t)dt = xi and / g(t) di = yt 

J i J i 



By Lemma 12.51 there exist closed intervals Ij = [aj,6j] C [0, m] where 1 < i < n such that for 
I = Ur=i[ a i) bi] it holds that 

1 



/(£) dt+ g(t) dt 

I [0,m]\J 



f(t)+g(t)dt. 



(2) 



[0,m] 



(i 



We may assume < ai < bi < 0,2 < b2 < ■ ■ ■ < a n < b n < m. For 1 < i < n let 

a'i := \ai + 1/2 J and 6< := |&i + V 2 J • 

Note that the a[, b[ are natural numbers such that < a^ < b^ < a' 2 < b' 2 < • • • < a' n < b' n < m. By 
the definition of / and g, for 1 < i < n it holds that 



J en 



(It 



+ 



fl(«) dt 



< z and 



/(*) dt 



b, 



+ 



g{t) dt 



<z, 



where \(vi, . . . ,V2 n )\ '■= ([f 1 1 , • • • , \v2n\) for i>i, . . . ,V2 n £ R. So if some Oj (resp., 6j) is replaced by 
a^ (resp., 60, then the left-hand side of ([2]) changes at most by z. Hence, for I' = Ur=i[ a i>^i] ^ 
holds that 

-2nz + ^ J f(t)+g(t)dt < ff(t)dt+f g(t)dt < 2nz+^ J f(t)+g(t)dt. (3) 

[0,m] I' [0,m]\I> [0,m] 

Let I" = Ur=iK> a- + 1, . . . , 6< - 1}. From Q and © we obtain 

m— 1 _, rn— 1 



-2nz + - ^(xi + yi) < ^ x» + ^J jfe < 2nz + - ^J (x, + 



■</; 



i=0 



is/" igJ" 



8=0 



D 



Next we state the integer variant of Lemma 12,61 

Corollary 2.7. Let n, rn > 1, xi,...,x m € Z 2n , and z G N 2n such that —z<Xi<z for all 
i. There exist natural numbers 1 < a\ < &i < 02 < 62 < • • • < a n < b n < m such that for 

1 = U"=l{ a i) Oj + 1, • • • , &i - 1}, 



-Anz < ^^ 



iei i<fi 



Xj < Anz. 



Proof. Let x\ := z + xi and y[ := z — x\. Thus x\,y\ € N 2n and x\,y\ < 2z. Lemma [231 applied 
to x\ and y[ provides natural numbers 1 < a\ < b\ < 02 < 62 < ■ ■ ■ < a n < b n < m such that for 

1 = ur=i{ a i' (H+i,...,bi-i}, 



-.in ..in 

-4nz + -^(x'i + y-) < J2 x 'i + J2 y i - 4nz + 2 53^ + y '^ 

iei i(£I 



i=\ 



i=l 



Therefore, 



2mz 1 v— v, v^ v~^ 2mz 1 ^— «.. 

-4nz + —— + - 2_^{xi - Xi) < mz + 2_^ x i~ /_^ x i ^ 4nz + — — + - } (xj - Xj). 

i=l iei i£I i=l 



a 



For the applications to maximum asymmetric traveling salesman and maximum weighted satisfi- 
ability we need the following variant of Lemma 12.61 While providing only a lower bound for the 
balanced sum, it estimates the rounding error more precisely. 

Lemma 2.8. Let n, m > 1 and xi, . . . , x m , yi, . . . , y m € N 2n . There exists an n' £ {0, . . . , n} and 
natural numbers I < a\ < b\ < 02 < b2 < • • • < a n > < b n i < m such that for I = Uieji n '}{ a ^ a * + 
l,...,bi}, 

y&i + 2/& 2 ^ 1- Vb n , + ^2 Xi + J2 Vi - 2 ^<( Xi + yi ^ 

Proof. Again let the indices of Xi and y, start at and define the componentwise integrable functions 
/ and g as in the proof of Lemma 12.61 So for i E {0, . . . , m — 1}, 



/■i+V 2 /•*+! 

/ f(t)dt = Xi and / g(t)dt = Vi . (4) 

Ji Ji+i/2 



By Lemma 12.51 there exists an n' € {0, ...,n} and closed intervals Ii = [ai,b{] C [0,m] where 
1 < % < n' such that for / = ljr=i[ a *> &*] ^ holds that 

Jf(t)dt+ J g(t)dt > I J f(t)+g(t)dt. (5) 

I [0,m]\I [0,m] 

Here we only need the inequality even though Lemma 12.51 states an equality. We may assume 

< a\ < b\ < a2 < 62 < • • • < a n i < b n > < in. (6) 

By the definition of / and g, the following holds for every % € {0, . . . , m — 1}: 

t€[*,i + i/2) => g(t) = (0,...,0) 

te[»' + VM + l) ==► /(t) = (0,...,0) 

Claim 2.9. We may assume that {aj,bj} % [i + 1/2, i + 1] and {bj,a,j+i} % [i,i + V 2 ] / or a ^ 
j € {1, ... , n'} and i G {0, . . . , m — 1}. 

Proof. If aj, 6j £ [i + 1/2, z + 1], then / is on [a,j, bj) and hence J" j f(t) dt = 0. Thus the left-hand 
side of ([5]) does not decrease if we remove the interval [cy, 6j] from /. Similarly, if fej, Oj+i S [i, £+1/2], 
then g is on [6j, Oj+i) and hence / 6 J+1 <?(£) (it = 0. Thus the left-hand side of © does not decrease 
if we replace the intervals [dj,bj] and [ey+i, 6j+i] by the interval [oj,6j + i]. Note that after these 
changes (which include a decrement of n'), ([6]) still holds. D 

Claim 2.10. We may assume that a\, . . . ,a n i G N and b\ + 1/2, . . . , b n i + X J2 € N. 

Proof. Assume ctj £ [i + 1/2, i + 1). By Claim 12.91 6j ^ [i + Y2, i + 1] and hence bj > i + 1. Since 
/ is on [i + 1/2, i + 1), the left-hand side of © does not decrease if we let a,j := i + 1. After this 
change, ([6]) still holds. 



Assume aj € [i,i + V 2 )- By Claim [2T9l bj-i $. [i,i + V 2 ] an d hence 6j_i < i (for j > 2). Since g is 
on [i, i + 1/2), the left-hand side of ([5]) does not decrease if we let dj := i. After this change, ([6]) 
still holds. 

Assume bj €[i + 1 /2, i + !)■ By Claim [279| a,- ^ [i + V2, i + 1] and hence a^ < i+ 1/2. Since / is on 
[« + V 2 ! * + 1)) the left-hand side of ([5]) does not decrease if we let bj := i + 1/2. After this change, 
© still holds. 

Assume 6j € [i,i + V 2 ) an d i < m. By Claim [2791 o J+ i ^ [«, % + 1/2] and hence Oj+i > i + !/2 (for 
j < n'). Since g is on [i, i + 1/2), the left-hand side of © does not decrease if we let bj := i + 1/2. 
After this change, (|6]) still holds. 

It remains to argue for the special case bj = m. By Claim 12.91 a j ^ [m — l /2, m] and hence 
dj < m — !/2. Since / is on [m — 1/2, m), the left-hand side of ([5]) does not decrease if we let 
bj := m — 1/2. After this change, ([6]) still holds. D 

If we split the integrals on the left-hand side of © according to I = \Jl = i[a,i,bj\, we obtain 

a\ n'-l I bi ai+1 \ m 

j g{t)dt + Y. [Jf( t ) dt + j 9{t)dt\ + J g(t)dt > \ J f(t) + g(t)dt. (7) 

i=1 V, b t I b n , [o,m] 

For i € {1, . . . ,n'} let a = h - 1/2. From Claim ETU] and © it follows that 

< ai < c\ < a% < C2 < ■ ■ ■ < a n > < c n i < m — 1. 

Together with @ we obtain: 



1 

g(t)dt = yo + yiH VVai-i 



b >. 

/ /(£) dt = x ai + x fli+ i H hie, 

Qi + 1 

/ g(t) dt = y Ci + y Ci+1 + ■■■+ y ai+1 -i 

bi 

in 

/ g(t) dt = y Cn , + y Cn , +1 H h y m _i 



In these sums, each index j € {ci, C2, • • • , c n >} appears exactly twice, once as Xj and once as yj. All 
remaining indices j G {0, . . . ,m — 1} \ {c\, C2, ■ ■ ■ , <v} appear exactly once, either as Xj or as y^-. 
Therefore, with I' = \J i£ n n n{ai,di + 1, . . . , q} the left-hand side of (|7|) is equal to 

y Cl +y C2 H Vy Cn , +^2xi + ^2yi. (8) 



Applying ^ to the right-hand side of ([7]) yields the desired inequality 



m—l 



t&r i£r j=o 

Corollary 2.11. Let n,m > 1 and xi,...,x m ,yi,...,y m ,z G N 2n such that y-i < z for all i. 

There exist n' < min(n, m) disjoint, nonempty intervals Ii,...,I n i C {1, . . . ,m} such that for 

J = JiU-UJ B » ; 

^ m 

n' ■ z + ^Xi + ^yi > -^{xi + yi). 

3 Applications to Mult i- Objective Optimization 

3.1 Preliminaries 

Consider some multi-objective maximization problem O that consists of a set of instances X, a 
set of solutions S(x) for each instance x G X, and a function w assigning a /c-dimensional weight 
w(x,s) G N fc to each solution s G <S(x) depending also on the instance i£l. If the instance x is 
clear from the context, we also write w(s) = w(x,s). The components of w are written as Wi for 
i G {1, 2, . . . , £;}. For weights a = (oi, . . . , a&), 6 = (&i, . . . , &&) G N fc we write a > 6 if a^ > bi for all 

*e{l,2,...,fc}. 

Let x G X. The Pareto set of x, the set of optimal solutions, is the set {s G S{x) \ Ss' G 
S{x) (w(x,s') > w(x,s) and w(x,s') ^ w(x,s))}. For solutions s,s' G S(x) and a < 1 we say s 
is a-approximated by s' if Wi(s') > a • Wi(s) for all i. We call a set of solutions a-approximate 
Pareto set of x if every solution s G S(x) (or equivalently, every solution from the Pareto set) is 
a-approximated by some s' contained in the set. 

We say that some algorithm is an a-approximation algorithm for O if it runs in polynomial time 
and returns an a-approximate Pareto set of x for all input instances i£l We call it randomized if 
it is allowed to fail with probability at most 1/2 over all of its executions. An algorithm is an FPTAS 
(fully polynomial-time approximation scheme) for a given optimization problem, if on input x and 
< e < 1 it computes a (1 — e)-approximate Pareto set of x in time polynomial in \x\ + Ye. If 
the algorithm is randomized it is called FPRAS (fully polynomial-time randomized approximation 
scheme) . 

3.2 fc-Objective Maximum Asymmetric Traveling Salesman Problem 

Definitions. Let k > 1. An N k -labeled directed graph is a tuple G = (V,E,w), where V is some 
finite set of vertices, E C V x V is a set of edges, and w : E — > N is a A;-dimensional weight 
function. We denote the i-th component of w by Wi and extend w to sets of edges by taking the 
sum over the weights of all edges in the set. A set of edges MC£is called matching in G if no two 
edges in M share a common vertex. A walk in G is an alternating sequence of vertices and edges 
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vo,e\,Vi, . . .e m ,v m , where V{ € V, ej € E, and ej = (vj-±,Vj) for all < i < m and 1 < j < m. 
If the sequence of vertices Vq, v\, . . . , v m does not contain any repetitions, the walk is called a path 
and if vo,v±, . . . ,v m -i does not contain any repetitions and v m = vq, it is called a cycle. A cycle 
in G is called Hamiltonian if it visits every vertex in G. For simplicity, we will interpret paths and 
cycles as sets of edges and can thus (using the above mentioned extension of w to sets of edges) 
write w(C) for the (multidimensional) weight of a Hamiltonian cycle C of G. 

Given some N fc -labeled directed graph as input, our goal is to find a maximum Hamiltonian cy- 
cle. We will also use the multi-objective version of the maximum matching problem. These two 
maximization problems are defined as follows: 

fc-Objective Maximum Asymmetric Traveling Salesman Problem 
(fc-MaxATSP) 

Instance: N -labeled directed complete graph (V, E, w) 
Solution: Hamiltonian cycle C 
Weight: w(C) 

fc-Objective Maximum Matching (fc-MM) 

Instance: N fc -labeled directed graph (V, E, w) 
Solution: Matching M 
Weight: w(M) 

Papadimitriou and Yannakakis [PYOO] give an FPRAS for fc-MM, which we will denote by 
£;-MM-Approx R and use as a black box in our algorithm. Since fc-MM-Approx R will be called multiple 
times, we assume that its success probability is amplified in a way such that the probability that 
all calls to the FPRAS succeed is at least l /2. 



High-Level Explanation of the Algorithm. We apply the balancing results to the multi- 
objective maximum asymmetric traveling salesman problem and obtain a short algorithm that 
provides a randomized i/^-approximation. This improves and simplifies the (1/2 — ^-approximation 
that was given by Manthey [Man09| . Essentially our algorithm contracts a small number of edges, 
then computes a maximum matching, adds the contracted edges to the matching, and extends the 
result in an arbitrary way to a Hamiltonian cycle. 

The argument for the correctness of the algorithm is as follows: Each Hamiltonian cycle H induces 
two perfect matchings (the edges with odd and the edges with even sequence number in the cycle). 
For each objective i, the weight of one of the matchings is at least 1 /2 ■ Wi(H). The balancing results 
assure the existence of a matching M such that for all objectives the inequality Wi{M) > 1/2 ■ wi{H) 
holds up to a small error. This matching can be approximated with the known FPRAS for multi- 
objective maximum matching. Moreover, by guessing and contracting a constant number of heavy 
edges in H our algorithm can compensate the errors caused by the balancing and by the FPRAS. 
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Figure 3: Contracting the edge e = (u, v) deletes all edges incident to v and sets the weights of every edge (u, x) to 
the weights of the edge («, x) for x £ V \ {u, v}. Any Hamiltonian cycle passes through some edge (u, x) and hence 
can be expanded to a Hamiltonian cycle through e by replacing (u, x) with the detour (u, v), (v, x). 

Contraction and Expansion of Paths. Suppose that for a given N fc -labeled complete directed 
graph G = (V, E, w) we wish to find some Hamiltonian cycle that contains a particular edge 
e = (u,v). This reduces to the problem of finding some Hamiltonian cycle in the N fc -labeled 
complete directed graph G" = (V ,E' ,w') where the edge e is contracted by combining the nodes 
u and v into a single node while retaining the ingoing edges of u and the outgoing edges of v. 
More formally, we remove v and all incident edges from G and set w'{u,x) = w(v,x) for every 
x £ V \ {u,v} (Figure [3]). Now suppose we find a Hamiltonian cycle C in G' . Then there exists 
some x such that (u,x) € C". Note that w'(u,x) = w(v,x). We replace the edge (u,x) in C" with 
the detour (u,v),(v,x) and obtain a Hamiltonian cycle C in G passing through e. Moreover, C 
preserves the weights of C in the sense that w(C) = w'{C') + w(e). 

The notion of edge contractions and expansions can easily be extended to sets of pairwise vertex 
disjoint paths (Figure H]), where each path is contracted edge-by-edge starting at the last edge of 
the path, and different paths can be contracted in an arbitrary order. We make this precise with 
the following definitions. 

Definition 3.1. Let G = (V,E,w) be some N -labeled complete directed graph, let (u,v) £ E, let 
P C E be a path uq, ei , m , C2, u-i, ■ ■ ■ , e r , u r , and let Q C E be a set of pairwise vertex disjoint paths 
Pi,P2,...,P r <^E. 

1. contract f u , v )(G) = (V \ {v},{e € E \ v is not incident to e},w'), where w'(x,y) = w(x,y) if 
x ^ u and w'(u, z) = w(v, z). 

2. contractp(G) = contract^ (contract,^ (. . . contract,,,, (G) . . . )) 

3. contractQ(G) = contract p ± (contract p 2 {. . . contract p r (G) . . . )) 

We sometimes identify a graph with its edge set and apply contract directly to sets of edges and not 
to graphs. In this case, we also interpret the value of contract as an edge set. 

Observe that the result of contracting several pairwise vertex disjoint paths does not depend on 
the order in which the paths are contracted. 

We define edge expansion in a similar manner. Note that Definition 13.21 becomes essential if G' is 
obtained from G by a contraction of some set Q of pairwise vertex disjoint paths in G. 
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Figure 4: Example for the contraction of the path {(u, v), (v, y)} in the graph G resulting in the graph G" 
and the following expansion of the tour {(u, x), (x, u)} in G" . The final tour in G includes the contracted 
path. 

Definition 3.2. Let G = (V,E,w) and G' = (V',E',w') be two ^H k -labeled complete directed 
graphs, let T C E' be a Hamiltonian cycle of G' , let (u, v) € E, let P C E be a path 
UQ,e\,ui,e2,U2, ■ ■ ■ ,e r ,u r and let Q C E be a set of pairwise vertex disjoint paths Pi, P2, . . . ,P r C 
E. 

1. expand/^ >v \ (T) = {(x, y) G T \ x ^ u} U {(n, v)} U {(v,x) \ (u,x) S T} 

2. expandp(T) = expand er (expand er _ 1 (. . . expand ei (T) . . . )) 

3. expandg(T) = expandp r (expandp rl (. . . expand Pl (T) . . . )) 

Again observe that the result of expanding several pairwise vertex disjoint paths does not depend 
on the order in which the paths are expanded. 

Proposition 3.3. Let G = (V,E,w) be some N -labeled complete directed graph, Q C E be a 
set of pairwise vertex disjoint paths, and G' = (V',E' ,w') = contractg(G). For any Hamiltonian 
cycle T' C E' of G' , the edges T = expandg(T') form a Hamiltonian cycle of G with w(T) = 
w'{T')+w{Q). " 

Approximation Algorithm. First we prove that the following algorithm computes a (1/2 — e)- 
approximation for /c-MaxATSP. Then Theorem 13.61 shows that a modification of the algorithm 
provides a 1 /2-approximation. 
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Algorithm: 2fc-MaxATSP-Approx R (V, E, w, e) 



Input : N 2fc -labeled complete directed graph G = (V, E, w) and even #V 
Output: set of Hamiltonian cycles of G 

1 foreach F C E with #F < 2k that is a set of vertex disjoint paths do 

2 G' := contractF(G); 

3 M : = 2A;-MM-Approx R (G',e); 

4 foreach M £ M do 

5 extend M in an arbitrary way to a Hamiltonian cycle T" in G'; 

6 output expandp(T'); 

7 end 

8 end 

Lemma 3.4. Let G = (V,E,w) be an N 2k -labeled complete directed graph with an even num- 
ber of vertices, e > 0, and T C E some Hamiltonian cycle in G. With probability at least 1 /2, 
2k-MaxATSP-Approx R (V,E,w,e) outputs a (1/2 — e)- approximation of T within time polynomial in 

\(V,E,w)\+ye. 

Proof Let G = (V, E, w) be an N 2fc dabeled complete directed graph with even m = jj=V, and let 
T be some arbitrary Hamiltonian cycle in G. 

Claim 3.5. There is a set F of vertex disjoint paths in T with j^F < 2k such that there is a 
matching M' in the graph (V',E' ,w') = contractp(G) with w'(M') > -^wiT) — w(F). 

Proof. We apply Lemma [2.8l to the sequence of edge weights of T. Having an even number of edges, 
we can write T sequentially as 

T = ui,ei,v±, fi,U2,e 2 ,v 2 , f2, ■ ■ ■ , u p ,e p , v p ,f p ,ui 

where Ui,V{ € V and ej, fa G T. 

Since w(ei),w(fi) € N 2k , Lemma 12.81 shows that there exist k! € {0, ...,k} and natural numbers 
1 < a\ < 61 < a 2 < b 2 ■ ■ ■ < ay < b k > < p such that for I = \Jie{i,...,k'}{ a ^ a i + !>•••> h}, 



1 P 

w(f bl )+w(f b2 ) + --- + w(f hl )+Y,He l ) + Y, w (^ ^ 2^ (w(6l)+U;(/l)) - (9) 

Let S = {fb t , fb 2 i- ■ ■ 1 fb k ,} U {ei I i G 1} U {fi \ i ^ I}. Observe that it is possible that adjacent 
edges are contained in S. Figure [5] gives an example. 

Observe that f or 1 < j < p and /o = f p the following holds. 

fr-i^eS ^^ 3i£ {l,...,k'},j = Oi (10) 

ejJjeS ^ 3ie{l,...,k'},j = h (11) 
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Figure 5: Some part of the cycle T, where SCT contains the depicted edges and is partially defined by hi = j, 
di+i —j + 2, bi+i - j + 3, and a i+2 = b i+2 = j + 4. 

Let F = {e ai , fb 1 , e 02 , /& 2 , . . . , e ay , /& , } and note that #F = 2k' . We argue that contracting F will 
transform any path in S into a single edge such that the resulting edge set is a matching: 

Suppose S contains some path P = {e r , f r , e r+ i, f r +i, ■ ■ ■ ,e s , f s }, where we assume P to be maximal 
(i.e., f r -i,e s+ i $l S). From (fl"0|) . (fTTj) . and ai < b\ < 0,2 < hi ■ ■ ■ < a^ < &&' we can draw the 
following conclusions: 

• e r , f r G P yields r = b% for some 1 < % < k' 

• f r , e r+ i, f r+ i G P yields r + 1 = a i+ i = 6 i+ i 

• / r+ i, e r+2 , fr+2 G P yields r + 2 = aj +2 = 6;+2 



• / s _i, e s ,f s eP yields s = a s _ r+i = b s - r +i 

Note that e r ^ F, since otherwise e& 4 G F, hence Oj = r, and by (flOl) . / r _i G S", which contradicts 
the maximality of P. Therefore, contracting {/&. , e ai+1 , /& i+1 ,...,/&.} C P transforms P into the 
single edge e r . A similar argumentation shows the same result for paths that start with some edge 
f r or end with some edge e s . Hence, contracting P transforms every path in S into a single edge, 
and M' = contract ir(S') is a matching in the graph (V' ,E' ,w') = contracts (G). We further obtain 

w'{M') = w{S) - w{F) 

>(f bl ) + w{f b2 ) + ■■■ + w(f hl ) + £>(*) + J>(/0 - w(F) 



w{ 



iei 



Hi 



1=1 



1 



w(T) - w(F) 



which proves the claim. 



□ 



We fix the iteration of 2/c-MaxATSP-Approx R where the algorithm chooses P as in the claim. By 
Claim [331 we know that there is a matching M' of G' = (V ,E' ,w') with 



w'(M') > -w(T) - w(F). 



(12) 
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Hence, with probability at least 1 /2 the set Ai contains some matching M of G' such that 

w'{M) > (l-e)w'(M') 

03/1 \ 

> (1 _ e) ^ w{T) _ w{F) j . 

We extend M to some Hamiltonian cycle T' of G' in an arbitrary way without losing weights. By 
Proposition 13.31 we can expand T' with F and obtain a Hamiltonian cycle T in G with 

w{f) = w'{T')+w{F) 
>w'{M) + w{F) 

>{l-e)(±w{T)-w(F)\+w{F) 

>(l-e) l - w (T) + ew(F) 

>{l-e)\w(T). 

Moreover, the running time of every operation of the algorithm, including the execution of the 
randomized maximum matching algorithm, and the number of iterations of the loops are polynomial 
in the length of the input and 1 /e, which completes the proof of the lemma. □ 

Theorem 3.6. /c-MaxATSP is randomized l /2-approximable. 

Proof. Let G = (V, E, w) be an N 2fc -labeled complete directed graph with even m = #V, and let 
T be some arbitrary Hamiltonian cycle in G. The proof can easily be extended to graphs with an 
odd number of vertices or objectives. 

For each 1 < i < 2k we choose an /; G T with w(fi) > —w(T) (i.e., fi is a heaviest edge of T with 
respect to component i). We let F C T be a smallest set with even cardinality containing all the 
fi. We get #F < 2k and 

w(F) > —w(T). (13) 

m 

F is a set of vertex disjoint paths in T and hence can be used to contract edges in G and T. Let 

G' = (y',E',w') = contract f(G) and T' = contract^ (T). Clearly, T' is a Hamiltonian cycle in G'. 

Moreover, Q has an even number of vertices. By Lemma 13.41 we can find in polynomial time a 

Hamiltonian cycle T in G' such that w'(T) > (1 — e)su? / (T / ), where e = 1 /m. Moreover, we can 

expand T with F to obtain a Hamiltonian cycle T in G such that 

w(f) = w'{f) + w(F) 

>{l-s)\w'{T') + w{F) 



l -w'(T') + l -w(F)^ + l -w(F) - £ -w'(T') 



= \w{T)+ l -{ W (F)-ew'(T')) 
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where (*) follows from 



w(F) - ew'{T') = w(F) w'{T') 

> w(F) w(T) 

m 

ED i i 

> — w(T) w(T) 

m m 

= 0. 

Note that, although we do not know the set F of heaviest edges in the Hamiltonian cycle, we can 
simply try all possible sets of heaviest edges, since the number of objectives 2k is constant. □ 



3.3 fc-Objective Maximum Weighted Satisfiability 

Definitions. We consider formulas over a finite set of propositional variables V. A literal is a 
propositional variable ueVor its negation v, a clause is a finite, nonempty set of literals, and a 
formula in conjunctive normal form (CNF, for short) is a finite set of clauses. A truth assignment 
is a mapping /: V — > {0, 1}. For some v G V, we say that I satisfies the literal v if I(v) = 1, and I 
satisfies the literal v if I(v) = 0. We further say that I satisfies the clause C and write 1(C) = 1 
if there is some literal I G C that is satisfied by /, and / satisfies a formula in CNF if / satisfies 
all of its clauses. For a set of clauses H and a variable v let H[v] = {C G H \ v E C} be the set of 
clauses that are satisfied if this variable is set to one, and analogously H[v] = {C G H \ v G C} be 
the set of clauses that are satisfied if this variable is set to zero. 

Given a formula in CNF and a /c-objective weight function that maps each clause to a ^-objective 
weight, our goal is to find truth assignments that maximize the sum of the weights of all satisfied 
clauses. 



fc-Objective Maximum Weighted Satisfiability (fc-MaxSAT) 

Instance: Formula H in CNF over a set of variables V, weight function w : H — >• N 

Solution: Truth assignment I: V — > {0, 1} 

Weight: Sum of the weights of all clauses satisfied by I, i. e., w(I) = Yl W (C) 

CeH 
I(C)=1 



High-Level Explanation of the Algorithm. We apply the balancing results to fc-MaxSAT. 
For a given formula H in CNF over the variables V, the strategy is as follows: Start with a list 
of the variables V and guess a partition of this list into 2k consecutive intervals. Assign 1 to the 
variables in every second interval and to the remaining variables. The balancing results assure 
the existence of a partition that yields an assignment whose weights are approximately one half of 
the total weights of H , up to an error induced by the variables at the boundaries of the partition. 
The error can be removed by first guessing a satisfying assignment for several influential variables 
V° of the formula. This results in a i/a- approximation for fc-MaxSAT. 
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Approximation Algorithm. We show that the following algorithm computes a !/2-approxima- 
tion for 2/c-MaxSAT. 



Algorithm: 2fc-MaxSAT-Approx(if, w) 
Input : Formula H in CNF over the variables V = {v±, . . . ,v m }, 2fe-objective weight 

function w: H — > N 2fc 
Output: Set of truth assignments I: V — > {0, 1} 

1 f oreach V° QV with #V° < (2k) 2 do 

2 let I(v) := for all veV°; 

3 G := {C G H | -3v G V° (v G C)}; 

4 y 1 := {t; g V \ V° | 2k ■ w{G[v}) £ w{H \ G)}; 

5 set I(v) := 1 for all v eV 1 ; 

6 V' :=V\{V°UV 1 ); 

7 f oreach oi, 6i, a2, &2, • • • ,afc, 6fc G {i | fj G V'} do 

8 f oreach u, G V do 

9 if 3j(aj < i < bj) then I{vi) := 1 else I{vi) := 

10 end 

11 output / 

12 end 

13 end 

Theorem 3.7. fc-MaxSAT is 1 /2-approximable. 

Proof. In the following, we assume without loss of generality that the number of objectives 2k is 
even. We show that the approximation is realized by the algorithm 2/c-MaxSAT-Approx. First note 
that this algorithm runs in polynomial time since k is constant. For the correctness, let (H, w) 
be the input where if is a formula over the variables V = {v±, . . . ,v m } and w: H —> N is the 
2/c-objective weight function. Let I : V — > {0,1} be an optimal truth assignment. We show that 
there is an iteration of the loops of 2fc-MaxSAT-Approx(.£f , w) that outputs a truth assignment I 
such that w(I) > w(I )/2. First we show that there is an iteration of the first loop that uses a 
suitable set V°. 

Claim 3.8. There is some set V° C {v G V \ I (v) = 0} with #V° < (2k) 2 such that for 
G = {C G H | ->3v G V° (v G C)} and anyveV\V° it holds that 

2k-w(G[v])^w(H\G) => I (v) = l. 

Proof. As a special case, if 4f{v G V \ I (v) = 0} < (2k) 2 , the assertion obviously holds for 
V° = {v G V | I (v) = 0}, since I (v) = 1 for all v G V \ V°. 

Otherwise, let V° = {u2kt+r \ r = 1,2, ... ,2k and t = 0,1, ... ,2k — 1}, where the U2kt+r £ V are 
defined inductively in the following way: 

(IB) H := H 
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(IS) 2kt + r - 1 -»• 2kt + r: 

— choose v G V \ {ui, . . . ,v,2kt+r-i} such that I (v) = and w r (H 2 kt+r-i\v]) is maximal 

- U 2 kt+r ■ = V 

— H2kt+r '■= H2kt+r-l \ H2kt+r-lM 

- a 2 kt+r :=w(52M+r-l[«]). 

We now show that the stated implication holds, so let v G V \ V° and j £{1,2,..., 2k} such that 
2k ■ Wj(G[v]) > Wj(H \ G). Because the union |Ji=i -^i— l[W] = H\G is disjoint, we get 



2k 2fc-l 2fc-l 

w( " 

r=l t=0 t=0 



and thus 

Wj ( G [v]) > ^ ^i. 

Hence, by a pigeonhole argument, there must be some t £ {0,1, ... ,2k — 1} such that Wj(G[v]) > 
(oc2kt+j)j- But since G C H 2 kt+j-i and thus even w^i^fct+j-i^]) > Wj(G[tJ]) > (a2fct+i)j, the only 
reason we did not choose v in the iteration 2kt + j (or even earlier) is that I (v) = 1. □ 

We choose the iteration of the algorithm where V° equals the set whose existence is guaranteed 
by Claim [3^81 Furthermore let G and V 1 be defined as in the algorithm and observe that by the 
claim it holds that I (v) = 1 for all v G V . Since I {v) = for all v G V°, the truth assignment I 
defined in the algorithm coincides with I on V° U V 1 . 

Let further V = V \ (V° U V 1 ) and G' = {C G G \ ^3v G V° (v G C) A ^3w G F 1 (v G C) A 3u G 
V (w G C V v G C)} be the set of clauses that are not yet satisfied by / but that could be satisfied 
by further extending /. 

Now we apply the balancing result. Let L' = V' U {tJ | v G V'}. For v i G V let 

w(C) ^ w(C) 



■I -i. 



#(CnL') and Vi- Is # (cni') 



and for v t G F° U F 1 let 

Xi = Vi = 0. 

It holds that 

^ a?i + yi = ^2 Xi + j/j = to(G'). 

Note that for all Vi G V , we have the bound yi < w{G'\vj\) < w(G[vi]) < j^w(H \ G) because of 
the definition of V' and V . Hence, for all V{ G V, 

yi<^w(H\G). 
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If we scale all values x% and j/j to natural numbers, then by Corollary I2.1H there exist k! < k 
disjoint, nonempty intervals J\, . . . , Jy C {1, . . . , m} such that for J = J\ U • • • U Jy it holds that 

X)x i + X)ift>^(G')-fc , ^u;(if\G?)>^KG / )- W (ff\G)). 

The algorithm tries all combinations of k (possibly empty) intervals J\ = [a\, b\}, . . . , J\~ = [o^, by]. 
In particular, it will test the combination of the k' nonempty intervals mentioned in Corollary |2.11i 
For / being the truth assignment generated in this iteration it holds that 

w({C G G' | 1(C) = 1}) > ^x J + ^ yi > \(w(G') -w(H\G)). (14) 

ieJ i<£j 

Furthermore, since I and I coincide on V \ V', we have 

w{{C £H\G'\ 1(C) = 1}) = w{{C £H\G'\ I {C) = 1}) (15) 

> w({C e H \ G | I (C) = 1}) 
= w({H\G}). (16) 

Thus we finally obtain 

w(I) = w({C €H\G' | 1(C) = 1}) + w({C G G 1 | 1(C) = 1}) 

m 

> w({C £H\G'\ 1(C) = 1}) + \(w(C) - w(H \ G)) 
^ w({C EH\G'\ I (C) = 1}) + \(w(G') - w(H \ G)) 

> \w({C eH\G'\ I (C) = 1}) + \w(G') 

> \w(I ). D 
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